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Abstract 

In this paper, we show that the L p -error of asymmetric linear approximation 
of the quadratic function Q(x) = Y^j=i x "] on simplices in M d of hxed volume is 
minimized on regular simplices. 

1 Introduction 

One important problem in geometry is to study the approximation (in a specified 
metric) of smooth convex bodies by various polytopes. For instance, questions 
of approximating convex bodies by inscribed or circumscribed polytopes, by 
polytopes with fixed number of vertices or faces, by polytopes of best approxi- 
mation, etc. have been studied in this direction. After some occasional results 
in the plane, the book of Fejes Toth [27J was the first to provide a large num- 
ber of problems, ideas and results on polytopal approximation in dimensions 
two and three, concentrating specifically on extremal properties of regular poly- 
topes. Many extensions have been made afterwards to higher dimensions, other 
metrics, etc. (see ,7, 8,[16l[T7] and references therein). 

On the other hand, the question on approximation of functions, defined on 
a polytope, by piecewise linear functions, generated with the help of triangula- 
tions, in L p -metrics is of great importance in approximation theory. The prob- 
lems of interpolation, best and best one-sided approximation of functions by 
linear splines (i.e., piecewise linear functions) have been considered. The ques- 
tion of finding optimal adaptive triangulations, i.e., triangulations that depend 
on the function being approximated (see, for instance, |13II14) V is of particular 
interest for applications. In order to construct adaptive triangulations many au- 
thors took into account the Hessian of the function (or curvature of its graph) 

(see oEUEHnanaizii)- 

Note that the construction of the best (in a specified sense) polytope for an 
individual convex body, or construction of the best triangulation for approxi- 
mation of a specific individual function, is an extremely difficult problem. The 
above-mentioned papers have dealt only with asymptotically optimal sequences 
of polytopes or asymptotically optimal sequences of triangulations. 

One possible method to construct an asymptotically optimal sequence of 
polytopes or triangulations is as follows. As the first step, we construct an in- 
termediate approximation of the convex body surface (or the function, respec- 



tively) by piecewise quadratic surface (function), and then solve the problem 
of approximating the piecewise quadratic surface (function) by piecewise linear 
ones. The latter, in turn (at least for d = 2), is equivalent to solving the fol- 
lowing optimization problem (we shall give its statement for approximation of 
functions) . 

Let a quadratic function Q and a simplex T in R d of unit volume be given. 
We shall consider the best L p -approximation of function Q by linear functions 
defined on T (or the best one-sided approximation, which coincide with the 
deviation of interpolant for positively definite Q), and the problem is to find a 
simplex T* , for which the corresponding error is minimal. (The known solutions 
of this problem are listed in Section 3.) 

Therefore, in a number of questions of geometry and approximation theory, 
it is important to find a simplex of fixed (unit) volume such that the error of the 
best approximation of a given quadratic function on this simplex in a specified 
metric or the best approximation with constraints (for instance, one-sided) is 
minimized. 

In approximation theory, there exists a tool to view both the problem of 
finding the best approximation without constraints and the problem of find- 
ing the best approximation with constraints "under one umbrella" . The latter 
can be viewed as the best approximation in the spaces with asymmetric norm 
or so-called (a, /^-approximation (see, for example, [3I21I2D]), when positive 
and negative parts of the difference between function and the approximant are 
"weighted" differently. Such type of approximations are of a separate interest, 
since they can be considered as the problems of approximation with non-strict 
constraints (sec below for more precise statements), when constraints are al- 
lowed to be violated, but the penalty for the violation is introduced into the 
error measure. We think that such a method could also be interesting for certain 
geometric problems. 

Therefore, the purpose of this paper is to prove the optimality of a reg- 
ular simplex in the problem of minimizing (over the simplices unit volume) 
the best (a, /3)-approximation in L p -metric of quadratic function 
Y2j=i x ] by linear functions. Note that with the help of linear transformations 
the solution of this problem allows us to obtain the solution of analogous opti- 
mization problems for an arbitrary positive definite quadratic 
form. 

The paper is organized as follows. Section 2 contains definitions, notations 
and rigorous statements of the problem we study and closely related problems. 
Previous results and the main result of the paper are given in Section 3. Section 4 
is devoted to the proof of the main result. 
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2 Notations, definitions and statements of the 
problems 

Let dsN and let M. d be the space of points x = (x\ , . . . , Xd) ■ Every point x G R d 
determines (and is determined by) the row- vector with coordinates (x%, . . . , Xd), 
and we shall reserve the notation x for such a vector. The Euclidean distance 
between points a, b G M d is defined as usually by 

1/2 



For a row- vector x, let x* be the column- vector transponent to x. For a square 
matrix J, we denote its transponent matrix by J*. 

For a measurable bounded set G C M. d , let L P (G), 1 < p < oo, be the space 
of measurable and integrable in the power p (essentially bounded if p — oo) and 
functions / : G — > K endowed with the usual norm 




\L P (G) 



/(/ G |/(x)|fdx)» ifl<p<co, 
1 esssup{|/(x)| : x G G} if p = oo. 



Let / G L p (G) and let a locally compact subset H C L p (G) be given. Set 
E(f;H)i, (Gi t° be the best approximation of the function / by H in the L p - 
norm, i.e., 

E(f;H) Lp(G y.= M{\\f-u\\ Lp{G) : ueH}. 

In addition, set 

E ± (f; H) Lp[G) := inf{||/ - u\\ Lp(G) : ±u < ±f, u £ (2-1) 

Quantity (|2.1[) is called the best approximation from above (E~ (/; H) Lp ^) 
or below (_E + (/; iJ) l p (g)) °f the function / by the subset if in the L p -norm. 
The quantities E~(f; H)l (g) an( i E + (f\ H) Lp ^ are called the best one-sided 
approximations. For a,/3 > and / G L P (G), let 

|/(x)| a>/9 :=a/+(x)+)8/_(x), 
where ff±(x) := max{±g(x); 0}. Define the asymmetric L p -norm as follows: 
II/IIW(G) := \n++f3.f-\\ LpiG) 

l(j G \f(x)Lpd*y ifl<p<co, 
v esssup{|/(x)| Qj ^ : x G G} i£ p = oo. 



Asymmetric norms in connection with various problems in approximation theory 
were considered in papers [TJ 21 [21] and books [TU [20l [22]. By 
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E(f; H) p -a t p denote the best (a, /^-approximation [3(20] °^ ^ ne function / by 
H in the Lp-norm, i.e., 

E(f; H) Lp;aAG) = inf {||/ - u|| W(G ) :ueH}. 

Note that for a = /3 = l, we have E{f\H) Lp , lAG) =E(f]H) Lp(G) . Babenko pQ 
proved that the following limit relations hold (see also [201 Theorem 1.4.10]): 

lim E(f;H) Lp . AAG) =E+(f;H) Lp{G) , 

lim E(f;H) Lp;aAG) =E-(f;H) Lp(G) . ^> 

This allows us to include the problem of the best unconstrained approximation 
and the problem of the best one-sided approximation into the family of problems 
of the same type, and consider them from a general point of view (for more on 
this motivation, see [H|3]). Because of the relation 

11/ - U K;1J3 = 11/ - W ~ 1)IK/ - P > !> 

the problem of the best (1, /3)-approximation can be considered as the problem 
of the best approximation with non-strict constraint / < u. This constraint is 
allowed to be violated, but the penalty 

G8"-l)||(/-u)_||£ 

for its violation is introduced into the error measure. In what follows, we shall 
allow the value +oo for a or /?, in that case identifying E(f;H) Lp a ^ G \ with 
the corresponding one-sided approximation. 
Let 

Sx(G) := {g(x) = ax* + c : a e M d , ceR,x£ G}. 

The space 5i(G) will be the main approximation set in this paper. Let also T — 
{t 1 , . . . , t d+1 } be the d-dimensional simplex with vertices t J , j = 1, . . . , d + 1. 
We shall consider the following optimization problem. 
Let Q(x) = xx*, and for T C R rf , set 

_ £(Q;Si(T)) w( t) 

<Tp; a ,p-Al ) ) 

| / | * 

where |T| stands for the d-dimensional volume of the simplex T . The purpose 
of this paper is to solve 

Problem 1. Find 

Op-u,l3;d '■= inf <J p] a,P;d{T) (2.3) 

and describe all simplices T, for which the inhmum in the right-hand part 
of (2.3) is achieved. 

A solution to Problem 1 will allow to solve the following related problems. 
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Problem 2. Find 



a p . 4 :=M ^'_; i+ r P (2-4) 



T 



and describe all simplices T, for which the infimum in the right-hand part 
of (2.4) is achieved. 



Problem 3. Find 



and describe all simplices T, for which the infimum in the right-hand part 
of (12.51) is achieved. 



3 History and the main result 

Note that the quantity E + (Q\S\(T))l p (T) coincides with the error of linear 
interpolation of the quadratic function Q(x) at the vertices of the simplex T, 
and the quantity E~(Q;S\(7'))l p (T) coincides with the error of tangential in- 
terpolation of Q(x) on the simplex T . In view of formulas (|2 . 2[) . we have 

a p;d= Km v P ;i,p-4 and a~ d = lim o> ;Q , M . 

The quantity a^. d has been considered for d — 2 in connection with the 
problem of finding the best triangulation Ajy consisting of N triangles of the 
set GcB 2 , provided that the L p -error of interpolation at the vertices of A^r 
of convex function / is minimized. 

The first attempt to find a^. d is due to D'Azevedo and Simpson [12], who 
computed cr+ . 2 . To the best of our knowledge, the progress on this problem can 
be outlined as follows: 

(1) d = 2, p = oo [12]; 

(2) d > 2, p = oo [26]; 

(3) d = 3,p = 2 0; 

(4) d = 2,p=l 0; 

(5) d = 2,p = 2 [IS]; 

(6) d>2,p eN [10]; 

(7) d = 2,p€ (0,oo) 0; 

(8) d > 2, p G (l,oo) HI]. 

Remark 3.1. Note that cr^ d for d = 2 and p G (1, oo) was independently found 
by Chen [11] and Babenko, ei al. [6]. 
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Remark 3.2. By L p -error in the case p <E (0, 1), we understand the following 
expression: 

E(f-,H) Lp(G y.= mi^Jjf(x)-uW\Vdxy : u e j . 

Remark 3.3. Infimum in Problem [3] is achieved only on regular simplices. 

To the best of our knowledge, Problem [2] was solved only in the case d = 2, 
p = 2 and a = (3 = 1 by Nadler [2"3"ll24j . In the next section, we shall give the 
solution of Problem 1 for all a, /3 > 0, 1 < p < oo and d £ N. 

Let To be a regular simplex of unit volume in M. d . The main result of our 
paper is the following: 

Theorem 3.1. Let a, (3 > 0, d £ N and 1 < p < oo. Then 

&V\u,P\d = CTp-a,P;d{To)- 

In view of (|2.2I) we obtain the following statements. 
Corollary 3.1. Let deN and 1 < p < oo. T/ien 

£(Q;<Si(To)) ip( To) 

°p5rf = ■ 

\To\ + " 

Corollary 3.2. Let d e N and 1 < p < oo. T/ien 



a ± 



p;d i'7-|i+ i 
|7o| " 

Recall that E + (Q;Si(To)) is the error of linear interpolation of 
function Q. 

The next section is devoted to the proof of this theorem. 

4 Proof of the main result 

Let a, (3 > be fixed throughout this section. Note that the value of the 
quantity E(Q; <Si(T))l j3 (j-) 1S independent of translations of the simplex T 
and its volume. For simplices T, T' C M d , we shall write T = T' if there exists 
a motion F of the space K d such that F(T) = T', and we shall write T ^ T' 
otherwise. 

The proof of the main theorem consists of two parts contained in the follow- 
ing two lemmas. 

Lemma 4.1. Let T be an arbitrary d-dimensional simplex of unit volume. Then 
there exists a constant C > 0, independent of T, such that 

a~p;a,0;d{7~) > C(diamT) 2 . 



G 



Lemma 4.2. IfT, T ^ 7o, is a simplex of unit volume in R d then there exists 
a simplex T* C M. d of unit volume such that 

0> ;Q ,/3;d(T) > Pp;a,P;d(T*). 

Indeed, in view of Lemma 4.1, there exists an optimal (i-dimensional simplex 
T' of unit volume such that a p a p-d — o- p . a p-d{T')- Then, Lemma [4.21 gives 
V = T . 

Proof of Lemma 4-1- Let Td '■= T = {t 1 , t 2 , . . . , t d+1 } be a simplex of unit 
volume. Assume that Ht 1 — "t 2 1 1 2 = diam7d- In addition, for j = 1, . . . , d — 1, 
let 7~ = {t\t 2 , . . . ,t 3+1 } be a simplex in W. 

First, let us consider the case 1 < p < 00. For j = 2, . . . , d, by /ij denoting 
the length of the height from the vertex t 3 of the simplex Tj to the simplex 
75-i. For every a = (ai, . . . , ad) £ K d , let a' = (ai, . . . , a^-i). Then, 



""o q B-dOd) = inf / 



ii.i 



= inf / / \u 2 + yy* - a d u - a'y* - c|£ * dydw 

h d 

h A / u \ 2p+(d-l) 



(I 



^ / 7T °%;*,p ; d-i(Td-i) du 



h d 

d P (T \ 

< a p;a,0;d-l\ ld —l)- 



2p + d 

Proceeding by induction on d, we verify that 

P ln -\ > hyhz . . .h d p 

t W;<A'«W ^ (2 J5 + 2)(2p + 3)---(2p + d) <T ^' 3 ' ll/lj 

> CP ■ h2k3 ■■■ hd ■ (diamT d ) 2p+1 = CP(diamT d ) 2p , 
a! 

where T is some positive constant independent of the simplex Td- 
Let us turn to the case p = 00. In this case, we obtain 

{T d )> (TO 

= inf sup \u 2 — ku ~ c\ a p > T(diam7d) i 

feeR,ceK ug [ 0idiamTd ] 



□ 
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Proof of Lemma 4-2- Let T = {w^w 2 ,!; 1 , . . . ,t d_1 } ^ To be a simplex of unit 
volume. Without loss of generality, we may assume that Hw 1 — t 1 )^ 7^ ||w 2 — 
t 1 )^. Clearly, we can always choose the coordinate system in R d so that the 
vertices of T have the following coordinates: 

w 1 = (-5,0,0,. ..,0), w 2 = (5,0,0,. ..,0), 

where 5 := ^Hw 1 — w 2 ||2, and the remaining vectors t^t 2 ,. 



following coordinates: 



t 1 =: 


( 




ai,i 














t 2 


( 




ai,2 


a-2,2 











t 3 


( 


b 3 


01,3 


0-2,3 


03,3 








t d-2 = . 


( 


bd-2 


ax,d-2 


d2,d-2 


a?>A-2 ■ 


• dd-2.d-2 





t d - 1 =: 


( 


bd-i 


ai,d-i 


&2,d-\ 


fl3,d-l • 


■ &d-2A-l 





have the 

), 
), 
), 



)• 

Note that in view of our assumption, b\ ^ 0. In addition, it can be easily seen 
that 



"3,3 

Let 



^ for all j = 1, . 



,d-l. 
b:= (&!,.. 



,&d-l) 



and 



A := 



Ai.i 01,2 
a 2 .2 




ai,3 

02,3 

03,3 




\ 
Since the matrix A is non-singular, set 

y = (vi, ■ ■ -,yd-i) 



&2,d-2 
a 3,d-2 



0>d-2,d-2 




ai,d-i \ 

fl2,d-l 
G3,d-1 



a.d-24- 



Let I be the identity matrix of size (d — 1) x (e? 



(4.1) 

1). In addition, set R := 



y y + I. It can be easily seen that the matrix R is positive definite. Therefore, in 
view of the Cholesky decomposition (a standard technique in numerical analysis, 
whose description can be found, for instance, in [15J, there exists an upper 
triangular matrix U = (uk,j)i<k,j<d-i such that 

R = U'U. 



Moreover, for every j = 1, . . . , d — 1, we have 



3,3 



Dj-i 



, where 



/l 



D := 1, D k := dct 



vl 



2/12/2 
2/12/3 



yiyk-i 

V 2/i2/fc 



2/12/2 
1 + 2/2 2 

2/22/3 



2/22/fc-i 
2/22/fe 



2/12/3 

2/22/3 

1 + 2/1 



2/32/fc-i 
2/32/fc 



2/i2/fc- 
2/22/fc- 
2/32/fc- 



2/i2/fc \ 
2/22/fc 
2/32/fc 



1 



2/fc 
2/fc-i2/fc 



1 2/fc-i2/fc 
1 + 2/fc/ 



(4.2) 
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for all k = 1, . . . , d — 1. Consequently, Ujj > for every j = 1, . . . , d — 1. 

Let Q = (qk,j)i<k.j<d-i be the diagonal matrix such that (j^j = Ujj, j = 
1, . . . , d — 1. We define 

U := Q X U. 

Then U is the unit upper triangular matrix. Therefore, 

R = U t Q 2 U. 

Set 

M := UA. 

Denote the elements of M by mk,j, i.e., M = {rrik,j)i<k,j<d-i- Note that for 
every j = 1, . . . , d — 1, it follows that rrijj — cijj. 

Let us now consider the simplex T — {w 1 , w 2 , t 1 , . . . , t d_1 }, whose vertices 
have the following coordinates: 



t 1 


■= ( 





ai,i 














t 2 


:= ( 





^1,2 


fl 2,2 











t 3 


:= ( 





"»i,3 




03,3 








ld-2 


:= ( 





mi,d-2 


m24-2 


m 3 ,d-2 ■ 


Cld-2,d-2 





ld-l 


:= ( 





mi,d_i 


m 2 ,d-i 


m 3id _i . 




Od-1, 



Obviously, the volumes of simplices T and T coincide. 

Now we shall construct the linear transformation 5 : R d — >• M d such that 
S'(r) = T. To that end, set 

h := blVT 1 e 

and define the linear transformation S with the help of the matrix 

/l h \ 


S= : (U)- 1 ' 

Note that dct S = 1. Then for every a = (ai, a 2 , . . . , a^) G K'' and c 6 1, it 
follows that 



L := Hxx* - ax* - c\\ Lp;aAr) = HvS^v* - aSv* - c|| (4.3) 
Let us consider the simplex T = {w 1 , w 2 , t 1 , . . . , t d_1 } such that 



t 1 


■= ( 


-h 


ai.i 














) 


t 2 


:= ( 


-b 2 


«1,2 


a2,2 











) 


t 3 


:= ( 


-b 3 


ffll,3 


a 2 ,3 


03,3 








) 


1d-2 


:= ( 


-b d -2 


Ol,d-2 


a2,d-2 


03,d-2 • 


• a,d-2,d-2 





) 


1d-l 


:= ( 


-bd-i 


Ol,d-l 




a3,d-l • 


■ ad-2,d-i 


Od-l,d-l 


) 
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It can be easily verified that the linear transformation S : M. d — > K d defined with 
the help of the matrix 

(I h \ 





s = 



(U)" 1 



/ 



transforms the simplex T into the simplex T. Therefore, for a := {—ax, 
a>2, • ■ ■ , ad), we obtain 



xx* - ax* - c| 



IvS^v* -aSv* - ell 



3("7T 



(4.4) 



Due to the symmetry of simplices T and T, we obtain L = L. Then 
from (|4.3p and (|4.4p . we derive that 



L = \{L + L)> 



Note that 



v(S t S + S t S)v t - (aS + aS)v* - 2c 



2D := S'S + S 1 ^ 



1 



h* [(U)" 1 ] 



V o 



1 

-h* 







h ^ 

(u)- 1 

/ 

/l h \ 





/I 





Vo 



[(u)- 1 ] 1 ; ; (u) 

... \ 



Since h = bM" 1 = bA _1 (U) _1 , we have 

h*h+ [(U)- 1 ]^!!)- 1 = [(UJ-^KA^^^bA- 1 +I](U)- 1 
= [(U)- 1 ]^^)" 1 = Q 2 . 

Therefore, 



D 



1 




Q 2 



/ 



and D is the diagonal matrix with elements 1,-Di, 757 ; - 
diagonal (numbers Dj, j = 1, . . . , d — 1, were defined in 



(4.5) 



LK, 



D d -2 



on the mam 
Let F be the 
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linear transformation defined with the help of the diagonal matrix F, having 
elements Df_ t , Df_ x ■ Df_ x ■ <J1=£ , Df_ x ■ ^ %^ on the main diagonal. 
It can be easily seen that det F = 1. Let T* = F~ l {T). Then, in view of (|4.5p . 

where g = — ^— (aS + aS)F and c' = — t— . Consequently, 

l 

0~p;a,/3;d(T) > D ^_ j Cr p;a ^ ;( J (T* ) . 

Let us show that the assumption Hw 1 — t H2 7^ ||w 2 — t || 2 yeilds 
Dd-i > 1. Indeed, since the matrix R is positive definite, it follows that 

Ai-i = det (y*y + 1) = (1 + A x ) • • • (1 + X d _ 1 ) > 1, (4.6) 

where Xj > 0, j — 1, . . . ,d — 1, are the eigenvalues of the matrix y*y (the 
existence of non-negative eigenvalues is guaranteed by positive semidefiniteness 
and symmetry of y'y). Inequality (|4.6j) becomes an equality if and only if Xj = 
for all j = 1, . . . , d— 1, or, that is the same, yj = for every j = 1, . . . , d— 1. In 
view of (|4.ip and non-singularity of the matrix A, this yields that bj = for all 
j = 1, . . . , d — 1. However, this contradicts the fact that b\ ^ 0, which follows 
from the assumption Hw 1 — t 1 ^ ^ ||w 2 — t 1 ^. □ 
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